Characteristics of plasma sheath in multi-component plasmas with three-ion species

The plasma sheath of a three ion species plasma is studied numerically, relying on the results of the experiment by Yip et al. (Phys. Plasmas 23:050703 (2016) to measure the positive ion velocities at the sheath edge. The positive ion species (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Ar^+$$\end{document}Ar+, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Kr^+,$$\end{document}Kr+, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Xe^+$$\end{document}Xe+) are assumed to be singly charged and to be characterized by the same temperature. It is shown that the sheath characteristics, viz. the particle number densities, the electrostatic potential and the space charge density profile in the sheath all depend on the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Kr^+$$\end{document}Kr+ concentration that is gradually added to the argon-xenon plasma as the third positive ion species. Also, the effect of ion-neutral collisions on the sheath properties is investigated numerically. Our results may be extended to a multi-ion plasma with more than two species of positive ions.

where v 0i and c si = (k B T e /m i ) 1/2 are the individual drift and sound velocity of each of the ion species, while n 0i and n 0e are the ion and electron densities at the sheath boundary, respectively. From (1), it is clear that unlike the case of a single species plasma where v 0i ≥ c s , the equality may be satisfied by speeds faster or slower than the individual sound speed corresponding to a (any) given ion species among the various ion species present in the plasma. Unlike an electron-ion plasma (i.e. with a single ion species), inequality (1) does not uniquely determine the flow speed of each ion species at the sheath edge if there is more than one species of ion. For instance, in a two ion species plasma, when the electron temperature far exceeds the ion temperature ( T e ≫ T i ), two possibilities exist for (1) to be satisfied 18 : either (i) all ions reach the sheath edge with a common velocity, called "the system sound speed" c sys = i n 0i k B T e /n 0e m i 1/2 , or (ii) each of the ion species has its own Bohm velocity at the sheath edge. Unlike the theoretical prediction 14 , experimental results for two-ion species plasmas did not confirm the latter solution but instead showed that the ion speed at the sheath boundary is closer to the system velocity, c sys , 4,[15][16][17][18][19][20][21][22][23][24] . In most of these experiments, the laser-induced fluorescence (LIF) technique was used to measure the velocity of each ion species at the sheath edge. More experimental evidence for the system sound speed solution has been provided by Lee et al. 15 . They studied the Bohm criterion in single and two ion species plasmas with LIF in Xe and Ar − Xe plasmas and showed that the argon and xenon ion velocities approach the ion sound speed of the system near the sheath boundary. Oksuz et al. 16 used electrostatic probes and ion-acoustic waves to measure the drift speed of ions in the plasma-boundary transition region. They measured the ion-acoustic wave speed at the sheath edge to be twice what it is in the bulk plasma in a two-ion-species plasma. None of these experiments provided a physical mechanism for these observations. On the other hand, theoretical and experimental investigations have shown that if the plasma contains multiple ion species, ion-ion two-stream instabilities can be expected to arise independently, in addition to ion-acoustic instabilities in the presheath under conditions of low neutral pressure and T e /T i ≫ 1 due to convection of one species onto another 17,19,[23][24][25] . Following these considerations, it appears that these flow-driven instabilities may play an important role in the determination of the ionic speeds at the edge of the plasma sheath. In 2009, Baalrud et al 17 showed that considering collisional friction enhanced by ion-ion streaming instabilities explained the discrepancy in the latter experiment. They showed that the collisional friction force causes the faster ion species to slow down and the slower species to accelerate. Considering this fact, they determined the solution of the plasma sheath criterion in two ion species plasmas. In fact, the theory suggested by Baalrud et al 17 indicated that the difference in velocities at the sheath edge depends on the relative concentrations of the two ions as follows: the difference in velocities is small, with both species approaching the bulk sound velocity, when the concentrations are comparable, and becomes large, with each species reaching its own Bohm velocity, when the relative concentration difference becomes significant. To test these findings, drift velocities of Ar and Xe ions were measured with the LIF technique 18 . The predictions were found to be in excellent agreement with the experimental data.
Although it might appear that, in the light of these experimental observations, the dynamics related to the Bohm criterion in multi-ion plasmas has been elucidated, a recent experiment by Yip et al. 26 showed that despite being one of the oldest problems in plasma physics, sheath formation and its criterion are still not fully understood. For the first time, they measured the drift velocities of two ion species at the sheath boundary of a three-ion component plasmas and found an unexpected feature regarding the ion speeds at the sheath edge, when there are three positive ion species: under most circumstances, the speed at which the ions enter into the sheath region matches neither the "system" sound speed nor their individual Bohm speed. Their experimental results demonstrated that if an additional third ion component was added to a two-ion plasma, the drift velocity (value) of each of the main two ion species at the sheath edge iterates from the system sound speed to their individual Bohm velocity as the concentration of the additional (third) ion species increases.
The aforementioned surprising findings of Yip et al. 26 work have motivated us to numerically investigate the sheath characteristics in an electropositive plasma with three ion species. We have undertaken a meticulous study of the spatial profile(s) of the (various) ion density, electrostatic potential and space charge density in a multicomponent plasma consisting of singly charged Ar + , Kr + and Xe + ions. Our results show that the presence of the third positive ion species increases the sheath potential but decreases the space charge and the sheath thickness. We combine recent results based on kinetic modeling with experimental results on three-component plasma to accurately compute the velocity of the ion species at the sheath edge which are the necessary initial conditions in our model. This allows us to investigate the effect of ion-neutral drag forces on the structure and to characterize the properties of the plasma sheath in a plasma containing three positive ions species, an opportunity offered for the first time, to the best of our knowledge.
This study is organized in four sections. Following this Introduction (Sect. 1), the basic equations of our model are laid out in Sect. 2. In Sect. 3, the model equations are solved numerically and the obtained results are analyzed. Our conclusions are eventually summarized in the final Sect. 4.

Model and basic equations
An unmagnetized electropositive plasma consisting of electrons and three types of singly charged positive ions with different masses and the same temperatures is considered. The electrons are assumed to be in thermal equilibrium, hence the electron density obeys the Boltzmann distribution 4,5,14 . Furthermore, the plasma is assumed to be in contact with a planar surface perpendicular to the x-axis and to have reached the steady state. Moreover, it is assumed that the physical parameters of the sheath vary only along the perpendicular direction to the surface (x-direction). www.nature.com/scientificreports/ Under the steady state condition and ignoring further ionization resulting in extra source/sink terms, the (stationary form of the) governing equations for each of the ion species will include the fluid continuity and momentum equations: where x is the distance from the sheath edge, n i , v i , m i , p i = n i k B T i and ν i = n n σ s,i v i are the density, velocity, mass, pressure and the effective ion-neutral collision frequency of the ith ion species ( i = 1, 2, 3 ), respectively, n n is the neutral gas density and σ s,i is the momentum transfer cross section for collisions between each ion species and neutrals. The coupled fluid equations are closed by Poisson's equation which reads: The quasi-neutrality condition imposed by the equilibrium requirement reads: where n 0i is the ith ion species density at the sheath edge ( x = 0).

It is convenient to introduce dimensionless variables
where c s1 is the sound velocity of the lightest ion species in this work and De = (k B T e /4πe 2 n 0e ) 1/2 . Using these variables, Eqs.
(2)-(5) can be rewritten in dimensionless form as follows: where N 01 = n 01 /n 0e , N 02 = n 02 /n 0e , N 03 = n 03 /n 0e , is a dimensionless collisionality parameter (that differs among the ion species). For our numerical analysis, we have adopted a standard 4th-order Runge-Kutta method with step-size h = 0.003 is used to solve Eqs. (6)-(9) with the boundary conditions u i (X = 0) = u 0i , N i (X = 0) = N 0i , φ(X = 0) = 0 and dφ/dX = 0 at the sheath edge and φ(X) = φ w on the wall. The differential equations are solved from the sheath edge (X = 0) towards the wall. To determine the wall position and hence the sheath thickness, we start from the sheath edge and move in space until the requirement φ(X) = φ w is fulfilled. It should be pointed out here that, following 26 , we have adopted the velocity values reported for the ions at the sheath edge. Therefore, flow-driven instabilities such as ion-ion two-stream instabilities were ignored in our computational model.

Results and discussion
We will investigate the sheath structure in an electropositive plasma consisting of electrons and singly charged Ar + , Kr + and Xe + ions with T e = 1.95 eV , T 1 = T 2 = T 3 = 0.03 eV , ϕ w = −90 eV and equal concentrations of argon and xenon ions (n 0Ar + = n 0Xe + ) , i.e. reproducing the experimental conditions of 26 . We assume that krypton is gradually added to the low pressure unmagnetized discharge into a mixture of argon and xenon. In the following, argon, krypton, and xenon ions are labeled as species 1, 2, and 3, respectively. Also, as mentioned in Sec. 3, there are three unspecified parameters α 1 , α 2 and α 3 in Eqs. (6)-(8) which depend on the collision cross section of the respective ion species. In the energy range of interest ( ∼ eV ), the ion-neutral cross section can be considered to be constant and is taken to be ∼ 5 × 10 −15 cm 2 for argon 1 . We have taken α 2 /α 1 ≈ 2.5 and and v 03 in our calculation. In addition, following the Yip et al. 26 work, we assume that krypton has its individual sound speed at the sheath edge ( v 02 = c s2 ), regardless of its concentration. As mentioned in 26 , this is a simple approximation based on the argument that the instability enhanced friction affects only the argon and xenon ions as the main cause of two-stream instability, i.e. not the krypton ions, which play no role in the occurrence of this instability. Arguably, this theoretical prediction should be tested experimentally.
Considering the above mentioned conditions and parameters, we examine the sheath structure in a multicomponent electropositive plasma containing three species of positive ions ( Ar + , Kr + and Xe + ) and electrons by studying the dynamical profile of the charged particle density, electric potential and space charge.  26 , we assume that the argon and xenon ions have the same concentration N 01 = N 03 . From Fig. 1, it is seen that the density distribution of the electrons as well as that of the net positive ions N + both decrease in space upon increasing the krypton concentration.
The evolution of the space charge σ = N + − N e inside the sheath region is depicted in Fig. 2 for three different values of the krypton concentration: N 02 = 0.15 , N 02 = 0.46 and N 02 = 0.61 . The other parameters are the same as in Fig. 1. Similar to Fig. 1, it is seen that an increase of the initial density of krypton ions causes the space charge in the sheath region to decrease. It is also observed that, with increasing N 02 , the density of electrons drops to zero more rapidly, whereas the positive ion densities drop slowly at the beginning. Therefore, the space  www.nature.com/scientificreports/ charge emerges at a peak value which indicates that in this region, more positive particles gather to shield the negative potential of the board. Figure 2 shows that this peak shifts towards the sheath edge with increasing N 02 . Moreover, it can be seen from this figure that the sheath thickness (values), the distance between the sheath edge (X = 0) and also the point where the condition φ = φ w is fulfilled, are all sensitive to the presence of krypton ions and in fact decrease as N 02 increases. Figure 3 shows the variation of the normalized electrostatic potential in the sheath region of Ar + − Kr + − Xe + plasma for two different concentrations of krypton ions and for the same parameters as in Fig. 1. It is observed that the normalized potential of the sheath decreases infinitesimally by increasing N 02 . This is consistent with the result of Fig. 2, which showed that the space charge decreases with increasing krypton concentration. The independence of the sheath width from the ion mass is another point seen in this figure, which is in agreement with the results of Fig. 2 as well as Child-Langmuir law 1 . Figure 4 shows the variation of the ion velocities throughout the sheath layer for different values of N 02 . The other parameters are the same as in Fig. 1. It is seen that an increase in Kr + concentration leads to an increase in the velocity of all ion species in the sheath region. Therefore, it is obvious from Figs. 1, 2, 3 and 4 that the krypton ion presence affects the profiles of the electrostatic potential, velocity and charged particle density distribution and the space charge of a three ion species plasma in a significant way.
In the light of the two main assumptions considered by Yip et al 26 , namely that the krypton velocity and the density of argon and xenon remain constant despite the increase in krypton concentration, which we have also adopted in our present work, but also in the absence of further experimental results, it would take more time to provide conclusive explanations for the changes observed as the krypton concentration increases; this task is therefore outside our scope for now.
We shall now investigate the effect of the ion-neutral drag force on the sheath dynamics. To do this, we will solve the basic equations of our model for different α 1 values. The effect of the drag force on the density distribution of the electrons and of the net positive ion species is shown in Fig. 5 for N 02 = 0.15 , v 01 = 1750 m/s , v 02 = 1500 m/s , v 03 = 1470 m/s , for different values of α 1 . This graph indicates that unlike the ion density distribution, the falloff of the density distribution of electrons through the sheath becomes faster as the drag force increases. In addition, similar to Refs. [27][28][29] , the decrease in sheath thickness due to the increase in α 1 is clearly visible in this figure. Moreover, it is seen from Fig. 5 that for α 1 < 0.001 the results are reduced to the collisionless case, within good approximation.
In Fig. 6 the dependence of the space charge density distribution on the ion-neutral drag force is demonstrated. From this plot, we see that an increase in the ion-neutral collision frequency leads to an increase in the space charge peak due to the reduction in positive ion velocity. In addition, by increasing the collisional effects (via increasing α 1 ), the charge accumulation position shifts towards the sheath edge. Figure 7 illustrates the effect of the collision force on the electrostatic potential. Similar to plasmas with two positive ion species 29 , it is seen that the electrostatic potential (the sheath thickness) increases (decreases) by increasing the ion-neutral drag force (increasing α 1 ). Taking into account the variation of the space charge with increasing α 1 , the behavior of φ in Fig. 7 is reasonable. We conclude that the dynamical properties of the sheath depend considerably on the ion-neutral drag force.
Finally, Fig. 8 shows the variation of the velocity of the ion species throughout the sheath layer with the drag force. From this figure, it is observed that the ion velocities decrease by increasing α 1 , which is in agreement with the results of previous works for plasmas containing either one or two ion species 10,11,27,29 . Moreover, it is seen that the lighter ion population has higher velocity towards the sheath wall than the heavier ion species.

Conclusion
Relying on the results of a recent experiment in attempt to evaluate the Bohm velocity of positive ions at the sheath edge of Ar + − Kr + − Xe + plasma, we have investigated the behavior of the electrostatic potential of the density distribution of the charged particles in the sheath region of an electropositive plasma consisting of electrons and three ion species consisting of argon, krypton and xenon positive ions. All types of ions were assumed to be singly charged and to be characterized by a finite temperature. Following the architecture of the experiment, we have gradually increased the krypton concentration in the argon-xenon plasma and we have    www.nature.com/scientificreports/ investigated the impact on the properties and on the structure of the sheath. Our results have shown that the presence of the krypton positive ion species affects the sheath structure and the density distribution profiles of the charged particles. It was found that an increase in Kr + concentration causes a decrease in the space charge density, which leads to a decrease in the sheath thickness and an increase in the sheath potential. It was also shown that the amplitude of the space charge density peak decreases and shifts towards the sheath edge as Kr + concentration increases. Moreover, we have shown that the ion-neutral drag force considerably affects the sheath potential and the density distribution of the charged particles. The space charge peak grows with increasing the collision force, and hence the sheath width decreases. In addition, collisional effects cause the amplitude of the space charge peak to increase and its position comes closer to the sheath edge. It was also shown that the velocity of the positive ion species in the sheath layer decreases by increasing the collision force.